Superradiance and instability of the charged Myers-Perry black holes 

in the Godel Universe 
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We consider scalar field perturbations of the asymptotically Godel 5-dimensional charged rotating 
black holes with two equal angular momenta. It is shown that the spectrum of proper oscillations 
of the perturbation includes superradiant unstable modes. The reason for the instability is the 
confining Dirichlet boundary condition at the asymptotically far region of the Godel Universe. The 
confining box makes superradiant modes extract rotational energy from the black hole and, after 
repeated reflections from the black hole, grow unboundedly. A similar instability takes place for 
rotating black holes in the asymptotically anti-de Sitter (AdS) space-time. 
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INTRODUCTION 



Stability of a black hole's space-time against small per- 
turbations is the basic requirement for its existence. Sta- 
bility acquires a special meaning in higher dimensional 
gravity, because there is no traditional uniqueness theo- 
rem for higher dimensional solutions with an event hori- 
zon. A number of "black" solutions with various topolo- 
gies have been found [1] and stability could be the cri- 
terium which could discard unphysical solutions 0, 
Through the gauge/gravity duality, classical instability 
of black holes can be interpreted as the onset of ther- 
modynamical phase transition in the dual field theory in 
some cases [4]. Therefore, stability of various black holes 
has been actively studied during past decade [5| . In most 
cases it is difficult to prove (in) stability of black holes 
analytically. Then, analysis of its proper oscillations, de- 
scribed by the called quasinormal modes, is used: If no 
unstable modes are found in the spectrum, the space-time 
is believed to be stable. 

Stability of various higher dimensional black holes have 
been proved in 0, Q in the context of string theory and 
higher dimensional gravity. Instability, once it happens, 
can be stipulated by different reasons. Various insta- 
bilities have been observed for higher dimensional black 
holes in |6[. A special kind of instability takes place for 
rotating black holes in the AdS space-time, because of the 
phenomenon called the superradiance |7[. The latter is 
the amplification of the reflected wave due to extraction 
of rotational energy from the black hole. If the perturba- 
tion propagate under the Dirichlet boundary condition at 
infinity, what occurs for asymptotically AdS black holes, 
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the superradiance naturally leads to instability. Superra- 
diance and the instability induced by it have been stud- 
ied in a great number of works (see for instance d, Q 
and references therein). It was shown [9[ that all (and 
only) superradiant modes are unstable for gravitational 
perturbation of the higher dimensional simply rotating 
asymptotically AdS black hole. 

In the previous paper [10] we considered the proper 
oscillation frequencies, called quasinormal modes, of the 
non-rotating 5-dimensional black hole immersed in the 
higher dimensional analogue of the Godel (rotating) Uni- 
verse [lll.[l2j|. The exact solution for the black hole metric 
was found by Gimon and Hashimoto in the supergravity 
12[ and was further investigated in a number of works 
13j. In [To| we have shown that the quasinormal spec- 
trum is similar to the spectrum of normal modes of the 
pure Godel space- time. The latter has a number of com- 
mon features with the spectrum of pure AdS space-time, 
where the scale of the universe's rotation j plays the role 
of the inverse anti-de Sitter radius. 

In the present research we shall consider quasinormal 
modes of a much more general space-time: 5-dimensional 
charged rotating black holes with both equal angular mo- 
menta in the Godel Universe. Our main result here is 
observation of the superradiant instability for such black 
holes. In comparison with the earlier work [10], we have 
two extra parameters, the black hole's rotational param- 
eters a and the black hole charge Q. We shall analyse the 
dependence of the quasinormal modes, and first of all of 
unstable superradiant modes, from these parameters and 
determine the parametric regions of instability. 

The paper is organized as follows: Sec II gives basic for- 
mulas for the 5-dimensional charged rotating black hole 
in the Godel Universe. Sec III is devoted to separation of 
variables in the test scalar field equation, Sec IV briefly 
describes numerical methods which we used for analysis 
of quasinormal modes. Finally, in Sec V we discuss the 
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obtained results and parametric regions of instability. 



II. CHARGED ROTATING ASYMPTOTICALLY 
GODEL BLACK HOLES 

The bosonic fields of the minimal (4+1)- supergravity 
theory consist of the metric and the one-form gauge field, 
which are governed by the following equations of motion 



Rfiv — 2 (^F jJia F v ^g^ v F ^ 

n fP^ v — _J_ fr aX^fiu pp pp . 



(1) 

(2) 



here, SaX-yfiv — \J det g^ v taX-y^u- 

In the Euler coordinates (£, r, #, 0), the solution for 
the 5-dimensional charged rotating asymptotically Godel 
black hole with two equal angular momenta a = a\ = <22 
is given by Q 



ds z = -f(r)dt 2 -q(r)raldt-h(r)r 2 (a 3 L ) 2 



dr 2 
v(r) 



where o\ 



+ — (dO z + di/; z + d<jf + 2 cos ###) , (3) 
/> + cosOdip, 



ty* fy*(J ry>0 



3jQa a 2 (M-Q) a 2 Q 2 



2r 4 4r 6 ' 

(4) 



h(r) = j 2 (r 2 + 2M+6Q)- 
f - 2M Q 2 

2M 8j(M + Q)(q + 2jM + 4jQ) 

= 1 - — + ~2 

2(M - Q)a 2 + Q 2 (l - 16ja - 8j 2 M - 24j 2 Q) 



Here M and Q are charge and mass of the black hole. 
When Q — and a — 0, the above solution is reduced 
to the Gimon-Hashimoto solution. When a = j = we 
have the 5-dimensional Reissner-Nordstrom solution. For 
Q = and j = we obtain the Myers-Perry black hole 
with two equal angular momenta. 



Perturbations of the scalar field in a curved background 
are governed by the Klein- Gordon equation 
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(5) 



Since the background metric has the Killing vectors d t , 
di/j, d(j,, one can choose the ansatz for the wave function 
as 

$(t, r, 0, 1>, 0) = e~ iujt + + im ^Y(e)R(r)r- 3 / 2 . 

(6) 

Substituting ((6]) into ((5]) and separating the variables, 
one can find that the angular part of the function satisfies 
the equation 



1 d . ri 2mn cos — m 2 — n 2 x . 

sin (9— + ; — t^~z + A Y(0) = 0, 



(?) 



sin OdO dO sin 2 

where A is the separation constant with the eigenvalues 

A = l{l + 1), i = max(|ra|, \n\) + i, i = 0, 1, 2 ... . 

Then, the equation for the radial part takes the wave- 
like form 

'| + QW)%) = 0, (8) 

where r* is the tortoise coordinate, which is defined as 

dr 



dr* 



v(r) 



(9) 



The effective potential can be written in terms of the 
coordinate r as follows [10j 



Q(r) = (1-4%)) (uj 



2mq(r) 



r(l -4/i(r)) 



(10) 



, x /4A 9 16ra 2 /i(r) 3v(r) + 6r*/(r) 



(l-4/i(r))r 2 



4r 2 



Now we are in position to perform analysis of the quasi- 
normal spectrum for the obtained wave equation. 



IV. BOUNDARY CONDITIONS AND 
NUMERICAL METHODS 



III. SEPARATION OF VARIABLES 

In order to derive the wave equation one can use the 
relation 

q 2 (r)+f(r)(l-4h(r))=v(r), 
which implies 
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r 6 sin 2 ( 



There are two positive solutions of the equation v(r) = 
0: the event horizon r + and the inner horizon r_ 

2r 2 ± = 2M-8ajM- 16j 2 M 2 -SajQ (11) 
- 48j 2 MQ - 32j 2 Q 2 ±Vp (r+ > r_ > 0), 

where 

P = (2M-8ajM-16j 2 M 2 -8ajQ-48j 2 MQ-32j 2 Q 2 ) 2 
-4(2a 2 M - 2a 2 Q + Q 2 - 16ajQ 2 - 8j 2 MQ 2 - 24j 2 Q 3 ). 
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FIG. 1: Real (blue) and imaginary (red) parts of Co as functions of a for m = 1, A = 2 (M = 1, Q = 0), j = 0.02 (left panel), 
j = 0.03 (middle panel), j = 0.04 (right panel). 



At the classical level, a wave cannot be emitted from 
the event horizon of a black hole, so that the standard 
boundary condition for a great variety of black holes is 
requirement of the purely ingoing wave at the event hori- 
zon, 



-iCor* y/l—4h(r*) 



2mq{r+) 



r* -oo, 



r + — 4r + /i(r+) 
With respect to the coordinate r (fT2j) reads 



uj-mVL+. (12) 



*oc(r-r + )-™ r^r+, 6 = ^ ~ 4 * (r+) . (13) 

t/(r+) 

At spatial infinity we have 

R(r -> oo) = C+^+(r) + C_^_(r), (14) 

where 



#±(r) = e ± *" r V ± ^1 



r 3 



with 



K = 2m 



4cjj 



(15) 
(16) 



+ 2juj(M + Q) (3 - 8ja - 16(M + 2Q)j 2 ) . 

Since the exponents e ±Ju;r2 have purely real index, they 
do not describe in-going or out-going waves. Therefore, 
we are unable to impose usual quasi-normal boundary 
conditions. However, we can use the analogy with AdS 
backgrounds and require Dirichlet boundary conditions 
at spatial infinity. This implies that C+ = for Re(juj) > 
or C- = for Re(juj) < 0. 

The non-trivial behavior of the functions ty± is ob- 
served when Re(uj) = 0. In this case both exponents 
have oscillatory behavior at spatial infinity. Thus, in or- 
der to impose the Dirichlet boundary conditions we must 



consider the factor r a± . When Re(uo) = one can find 
that 



Re(a±) 



1 



±2ra, 



implying that one of ty± is convergent and the other one 
is divergent as r — » oo. The only exception is m = 0, 
when both and are convergent. However, in this 
case the function norm is divergent, so we conclude that 
such modes do not exist [loj . 

Now, we shall briefly relate the two standard meth- 
ods used for numerical search of quasinormal modes: the 
shooting method and Frobenius method. These two al- 
ternative methods were used in order to guarantee the 
validity of the obtained results. More detailed discussion 
of both methods can be found in [3[ . 

Shooting method. Since at the horizon we require 
the purely in-going wave (fT3]h it is convenient to define 
a new function in such a way, 



y(r) = 



ioob 



x R(r), 



(17) 



that it becomes regular at the event horizon. 

We fix the wave- function norm so that 2/(r+) = 1. 
Substituting (fT7|) into (|8j) and expanding the wave equa- 
tion near the horizon, we find that ?/ / (r + ). This gives us 
boundary condition at the horizon for any fixed u, which 
we use for the numerical integration of the equation ([5]). 
At large distance we compare the result of our numerical 
integration with the large-distance asymptotic series ex- 
pansion (fT4|) and find the coefficients C± by using the fit- 
ting procedure. Then, quasinormal modes can be found 
by minimizing C+(co>) or C-(uj) depending on the sign of 
Re{uo). In order to check convergence of the procedure 
we check that obtained frequencies do not change within 
specified accuracy, if we increase precision of floating- 
point operations or the distance at which we use the fit. 

Frobenius method. This method allows us to find 
QNMs by solving numerically an equation with contin- 
ued fractions, which takes much less computer time for 
finding quasinormal modes. One can rewrite the wave- 
function as 



R(r) 



-iCob 



x z(r), (18) 
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FIG. 2: Region of instability (blue) for (neutral) Myers- Perry- 
Godel black holes. 
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FIG. 4: Real part of uj as functions of a for m = 1, A = 2 (blue, 
bottom), m = 2, A = 6 (green), m — 3, A = 12 (orange), 
m = 4, A = 20 (red), m = 5, A = 30 (magenta, top) (M = 1, 
Q = 0,j = 0.03). 
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FIG. 3: Region of instability (blue) for charged rotating 
asymptotically Godel black holes for j — 0.02 (M=l). 



where z(r) must be regular at spatial infinity and the 
event horizon, once uj is the quasinormal frequency. We 
choose "-" sign for Re(ujj) > 0, and sign for 

Re(ujj) < 0. 

The function z(r) can be expanded into a series near 
the horizon 

00 / r 2 _ r 2 \ i 

*M = £M;^ )■ ( 19 ) 

2=0 V - / 

After we substitute ([T8]) and ([T9]) into equation (j5J), we 
find the five-term recurrence relation for the coefficients 
di. Using Gaussian eliminations we numerically reduce 
the five-term recurrence relation to the three-term recur- 
rence relation [3j and solve the equation with the infinite 
continued fraction with respect to uj [l6[ . 

V. SUPERRADIANT INSTABILITY 

Similarly to the asymptotically AdS rotating black 
holes, the superradiance occurs when 

Re(uj) < mQ + . (20) 



On Fig. [T]in the regime of small j, we have shown that 
only the superradiant modes are unstable, while non- 
superradiant ones are stable. One can also see on Fig. [1] 
that the larger j corresponds to the higher superradiant 
instability growth rate. The highest growth rate which 
we have observed is about ~ 10 -6 M -1 / 2 . 

For j = positive and negative a produce the same 
quasinormal spectra (under the change m — » —m) be- 
cause of the identical picture of "left" and "right" rota- 
tions. On Figs. [2] and [3] one can see that this symmetry 
between left and right rotations is broken due to the ro- 
tation of the universe. Now, negative a allows for larger 
parametric region of stability than positive one. In addi- 
tion, increasing of the Universe rotation "stabilize" the 
system, so that already for jy/M ~ 0.075 there is no 
superradiance for any a (Fig. [2j). 

It is well known that the Reissner-Nordstrom met- 
ric depends only on square of the black hole charge 
Q. Therefore, the Reissner-Nordstrom space-time is the 
same for positive and negative charges. Unlike, the RN 
case, the Gimone-Hashimoto metric and its rotating gen- 
eralization are not "symmetric" respectively the change 
Q — >• —Q. The dependence of the QNMs on the charge Q 
is also not "symmetric" for opposite signs of the charge. 
Larger values of Q allows for bigger parametric "island" 
of stability, so that the near extremal positively charged 
black hole look like having no superradiance at any non- 
extremal rotation parameter a (Fig. [3]). In addition, 
near extremal negative values of Q evidently correspond 
to stability. In the regime of large Universe scale j, the 
quasinormal frequencies approach the normal modes of 
the pure Godel space-time [10|, LLZl similarly to the spec- 
trum of AdS balck holes in the limit of vanishing AdS 
radius [lEj . 

As the superradiance condition ([20]) is almost equidis- 
tant respectively m for small j, the region of instability 
practically coincides for all m because Re(uj) is equidis- 
tant with respect to m for small j. Therefore, as can be 
seen on Fig. HJ when determining the threshold of insta- 
bility, one can be limited by lower m within acceptable 
accuracy. 



5 




FIG. 5: Real (left panel) and imaginary (right panel) parts of the dominant quasinormal mode as a function of a and Q for 
j 0.02 (ra = 1, A = 2, M = 1) 




FIG. 6: Real (left panel) and imaginary (right panel) parts of the dominant quasinormal mode as a function of a and Q for 
j 0.07 (m = 1, A = 2, M = 1) 



Quasinormal modes as function of the black hole 
charge Q and the rotation parameter a are shown on 
Figs. [6j[7J[8j[9l When increasing the black hole charge Q 
the real oscillation frequency of the quasinormal mode 
uo monotonically increases (see Figs. [5j [6]), while the 
damping rate usually decreases, making QN modes of 
the charged black hole longer lived and having higher os- 
cillation frequency. This means that charged black hole 
is a better oscillator (i.e. has the larger quality factor, 
which is proportional to \Re(u) / Im(uS)\) than the neutral 
one. A similar dependence of the QNMs on the param- 
eter of the black hole's rotation a takes place, once the 
other parameters (M, Q, j) are fixed. Therefore, in the 
region of stability, the charged rotating asymptotically 
Godel black hole is considerably better oscillator than 
the neutral non-rotating one. 



VI. CONCLUSIONS 

In the present work, we have considered scalar field 
perturbations of the asymptotically Godel 5-dimensional 
Myers-Perry black holes with two equal angular mo- 
menta. Due to the Dirichlet boundary conditions of the 



wave equation which are required for the asymptotically 
Godel space-time, the quasinormal spectrum has similar 
features to the 5-dimensional Myers-Perry black holes in 
the anti-de Sitter space-time. A superradiant waves for 
such a rotating black hole, when immersed in an effective 
confining box (i.e. under Dirichlet boundary conditions) 
leads to the superradiant instability for some values of the 
black hole parameters. A similar kind of instability we 
have observed for the Myers-Perry- Godel black holes. It 
is evident to us that the same kind of instability should 
occur also for gravitational perturbations, i.e. Myers- 
Perry-Godel black holes must be unstable. Further, we 
have accurately determined the parametric region of in- 
stability for various values of the Universe rotation, black 
hole rotation and the black hole charge. The found in- 
stability has a relatively small growth rate. 

In our opinion, the instability found here for the partic- 
ular black hole solution (the five-dimensional black hole 
in the presence of the gauge field) is essentially generic 
phenomena which is induced only by the two factors: 

1. the rotation of the black hole, which leads to su- 
perradiance and 

2. the Dirichlet boundary condition in the asymptotic 
region of the rotating Universe. 



FIG. 7: Real (left panel) and imaginary (right panel) parts of the dominant quasinormal mode as a function of a and Q for 
j = 0.07 (m = — 1, X = 2, M = 1) 




Therefore, we expect similar instability for any other 
types of rotating black holes immersed in a rotating 
Godel-like Universe. 

Question which was not completely solved here is the 
stability of the extremely (positively) charged rotating 
black holes: Although numerical data indicates that the 
range of values of a for which instability takes place 
shrink to zero when approaching the extremal positive 
charge Q, there is no analytical proof for this observa- 
tion. 

Another interesting question is the dynamical evolu- 
tion of perturbation for the above case, which is quite 
non-trivial problem due to the existence of the closed 
time-like curves far from the black hole. There, the coor- 
dinate t becomes spacelike and the azimuthal coordinate 
(j> is timelike. This problem arises if one considers evolu- 
tion of the perturbation in time domain. In the frequency 
domain we are free from this kind of problems and, there- 
fore, we have considered here only self-consistent solu- 



tions. 
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